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ABSTRACT 

Axisymmetric pulsations of rotating neutron stars can be excited in several 
scenarios, such as core-collapse, crust and core-quakes and binary mergers 
and could become detectable either in gravitational waves or high-energy 
radiation. Here, we present a comprehensive study of all low-order axisym- 
metric modes of uniformly and rapidly rotating relativistic stars. Initial 
stationary configurations are appropriately perturbed and are numerically 
evolved using an axisymmetric, nonlinear relativistic hydrodynamics code, 
assuming time-independence of the gravitational field (Cowling approxima- 
tion). The simulations are performed using a high-resolution shock-capturing 
finite-difference scheme accurate enough to maintain the initial rotation law 
for a large number of rotational periods, even for stars at the mass-shedding 
limit. Through Fourier transforms of the time evolution of selected fluid vari- 
ables, we compute the frequencies of quasi-radial and non-radial modes with 
spherical harmonic indices I = 0, 1, 2 and 3, for a sequence of rotating stars 
from the non- rotating limit to the mass-shedding limit. The frequencies of 
the axisymmetric modes are affected significantly by rotation only when the 
rotation rate exceeds about 50% of the maximum allowed. As expected, at 
large rotation rates, apparent mode crossings between different modes ap- 
pear. In addition to the above modes, several axisymmetric inertial modes 
are also excited in our numerical evolutions. 

Key words: Hydrodynamics — relativity — methods: numerical — stars: 
neutron — stars: oscillations — stars: rotation 



1 INTRODUCTION 

The pulsations of rotating neutron stars are expected 
to be a source of detectable gravitational waves. Ad- 
ditionally, their excitation could become detectable in 
the emission of high-energy radiation. In particular, 
axisymmetric oscillations can be excited in a number 
of different astrophysical scenarios, namely: a) after a 
core-collapse leading to a supernova explosion (see e.g. 
Monchmeyer et al 1991; Zwerger & Muller 1997), b) 
during starquakes induced by the secular spin-down of 
a pulsar, c) after a large thermonuclear explosion in the 
crust of an accreting neutron star, d) during a core- 
quake due to a large phase-transition to, for example, 
strange quark matter (Cheng & Dai 1998), and e) in 



the delayed collapse of the merged object in a binary 
neutron star merger (Ruffert, Janka & Schafer 1996; 
Shibata & Uryu 2000). The observational detection of 
such pulsations will yield valuable information about 
the equation of state of relativistic stars (see Kokko- 
tas, Apostolatos & Andersson, 2000, see also Kokkotas 
and Schmidt 1999, for a recent review on oscillations 
of relativistic stars) . 

Numerical simulations of some of these scenarios 
are available and provide very detailed information of 
the dynamics of the neutron star pulsations. In par- 
ticular, the axisymmetric core-collapse simulations of 
Monchmeyer et al (1991) and Zwerger & Muller (1997), 
revealed that, after the collapse and bounce of an iron 
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core, the unshocked inner core (the proto-neutron star) 
oscillates with various volume (radial and quasi-radial) 
and surface modes. The amplitude and frequency of 
these fluid modes (/- and p-modes) was found to de- 
pend on the kinetic energy of the inner core at bounce, 
the stiffness of the equation of state (EOS), and the 
central and average density of the inner core. These 
authors found that the amplitude of the post-bounce 
oscillations is small for spherical models, being strongly 
damped through the emission of asymmetric pressure 
waves, in time scales of the order of 1ms. However, for 
rotating cores which bounce due to centrifugal forces 
at subnuclear densities, much larger amplitudes are 
achieved (as large as 10 times the central density) and 
the damping time scale becomes comparable to the os- 
cillation time scale (S> 1ms). Recently, Dimmelmeier, 
Font & Miiller (2000) have developed a code to study 
axisymmetric core-collapse in general relativity using 
the conformally flat metric approach (Wilson, Mathews 
& Marronetti 1996). This code is currently being ap- 
plied to collapse some of the initial models of Zwerger & 
Miiller (1997), to analyze the gravitational waves emit- 
ted in the process. Excitation of axisymmetric modes 
have already been observed in such relativistic core- 
collapse simulations. 

Of all axisymmetric modes, the quasi-radial modes 
of slowly-rotating relativistic stars were first studied by 
Hartle & Friedman (1975) and, more recently, by Datta 
ct al. (1998). In rapid rotation, quasi-radial modes 
of relativistic stars have been studied by Yoshida & 
Eriguchi (2000) in the Cowling approximation (Mc- 
Dermott, Van Horn & Scholl 1983), i.e., by neglecting 
the perturbations in the gravitational field (see Ster- 
gioulas 1998 for a recent review on the equilibrium 
structure and oscillations of rapidly rotating stars in 
general relativity). For Newtonian stars, axisymmetric 
modes have been extensively studied by Clement (1981, 
1984, 1986). From the above studies, it has become 
apparent that rotation weakly modifies the oscillation 
frequencies for low-order modes, but introduces appar- 
ent crossings between higher-order modes for rapidly 
rotating models. In addition, in Clement (1981) it is 
claimed that the axisymmetric quadrupole /-mode lies 
on a continuous branch for rapidly rotating Newtonian 
stars. 

In this paper we compute all low-order I = 0, 1, 2 
and 3 axisymmetric modes for rapidly rotating stars in 
general relativity, in the Cowling approximation. For 
this purpose, we use a 2-D nonlinear hydrodynamics 
code, whose accuracy has been extensively tested in 
Font, Stergioulas & Kokkotas (2000) (hereafter FSK; 
see also Stergioulas, Font & Kokkotas 1999). This code 
is based on high-resolution shock-capturing (HRSC) 
finite-difference schemes for the numerical integration 
of the general relativistic hydrodynamic equations (see 
Font 2000 for a recent review) . We note in passing that 
the 3-D version of the numerical methods employed 



here has been recently applied by Stergioulas & Font 
(2000) in the study of the large-amplitude, nonlinear 
evolutions of r-modes in rotating relativistic stars. In 
our present study of axisymmetric modes, we focus 
on a sequence of equilibrium models with a polytropic 
(N — 1.0) equation of state and uniform rotation. For 
the excitation of the various oscillation modes, low- 
amplitude perturbations (using appropriate trial eigen- 
functions) are added to the initial equilibrium models. 
The Cowling approximation allows us to evolve rela- 
tivistic matter for a much longer time than presently 
available coupled spacetime plus hydrodynamical evo- 
lution codes (Alcubierre et al. 2000; Font et al. 2000; 
Shibata, Baumgarte & Shapiro 2000; Shibata & Uryu 
2000). This is particularly evident when hydrodynam- 
ically evolving rotating stars. Nevertheless, since pul- 
sations of neutron stars are mainly a hydrodynamical 
process, the approximation of a time-independent grav- 
itational field still allows for qualitative conclusions to 
be drawn, when studying the evolution of perturbed 
rotating neutron stars. In addition, our present results 
will serve as test-beds for 3-D general-relativistic evo- 
lution codes. 

The paper is organized as follows: In Section 2 we 
describe the setup of the problem, by briefly presenting 
some details of the initial equilibrium stellar configu- 
rations and the main features of the hydrodynamical 
code. In Section 3 we explain the procedure by which 
the initial equilibrium models are perturbed. Section 4 
presents the main results of our simulations, including 
the frequencies of all low-order axisymmetric modes for 
our sample of initial models. The paper ends with Sec- 
tion 5 where a summary is presented, together with an 
outlook of possible future directions of this investiga- 
tion. 



2 PROBLEM SETUP 

Our initial models are numerical solutions of the exact 
equations describing rapidly rotating relativistic stars, 
having uniform angular velocity f2. We assume a per- 
fect fluid, zero-temperature equation of state (EOS), 
for which the energy density is a function of pressure 
only. The following relativistic generalization of the 
Newtonian polytropic EOS is chosen: 

V = K P \^\ (1) 
e = po + Np, (2) 

where p is the pressure, e is the energy density, pa is 
the rest-mass density, K is the polytropic constant and 
N is the polytropic exponent. The initial equilibrium 
models are computed using a numerical code devel- 
oped by Stergioulas & Friedman (1995). Our represen- 
tative neutron star models are characterized by N = 1, 
K = 100 and central density p c = 1.28 x 10~ 3 , in 
units in which c — G = M Q = 1. We compute 12 
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Table 1. Equilibrium properties of the initial 
models, as described by a polytropic EOS, p = 
Kp^ 1 ^ , where N = 1, K = 100 and with central 
rest-mass density p c = 1.28 X 10~ 3 (in units with 
c = G = Mq = 1). The entries in the table are 
as follows: fi is the angular velocity of the star, M 
and Mo are the gravitational and rest mass, T/W 
is the ratio of rotational to gravitational binding 
energy and R is the equatorial circumferential ra- 
dius. 



n 


M 


M 


T/W 


Ft 


(io 4 ^- 1 ) 


(M ) 


(M Q ) 


(xl0~ 2 ) 


(km) 


0.0 


1.400 


1.506 


0.0 


14.15 


0.218 


1.432 


1.541 


1.200 


14.51 


0.306 


1.466 


1.579 


2.438 


14.92 


0.371 


1.503 


1.619 


3.701 


15.38 


0.399 


1.523 


1.641 


4.339 


15.63 


0.423 


1.543 


1.663 


4.976 


15.91 


0.445 


1.564 


1.686 


5.609 


16.21 


0.465 


1.585 


1.709 


6.232 


16.52 


0.482 


1.607 


1.733 


6.839 


16.87 


0.498 


1.627 


1.756 


7.419 


17.25 


0.511 


1.647 


1.778 


7.959 


17.68 


0.522 


1.666 


1.798 


8.439 


18.15 



different initial models by varying the polar to equa- 
torial circumferential radius from 1 (non-rotating star) 
to 0.65 (near the mass-shedding limit), as listed in Ta- 
ble 1. The angular velocity at the mass-shedding limit 
is Qk = 0.8094 x 10 4 s _1 for this sequence of rotating 
relativistic stars of same central density. In order to be 
able to study stellar pulsations the initial model is sup- 
plemented by a uniform, non-rotating "atmosphere" of 
very low density, typically 10~ 6 or less times the central 
density of the star (see related discussion in FSK) . 

The initial data are subsequently evolved in time 
with a hydrodynamics code. The (axisymmetric) hy- 
drodynamic equations are written as a first-order flux- 
conservative system which expresses the conservation 
laws of mass, momentum and energy. The specific form 
of the equations was presented in FSK and we will 
not repeat it here. These equations are solved using a 
HRSC finite-difference scheme (Font 2000). A compre- 
hensive description of the specific numerical techniques 
we employ was previously reported in FSK. Therefore, 
we only mention here that the code makes use of the 
third-order piecewise parabolic method, PPM (Colella 
& Woodward 1984), for the cell-reconstruction proce- 
dure, together with Marquina's flux-formula (Donat et 
al 1998) to compute the numerical fluxes. The PPM re- 
construction scheme was shown to be accurate enough 
for maintaining the initial rotation laws during many 
rotational periods. 



The hydrodynamic equations are implemented in 
the code using spherical polar coordinates (r, 0, 4>) and 
assuming axisymmetry, i.e., all derivatives with respect 
to the 4> coordinate vanish. The radial computational 
domain extends to 1.2 times the radius of the star (the 
20% additional zones are used for the atmosphere). In 
the polar direction, the selected domain depends on 
the spherical harmonic index of the pulsation modes: 
for even I modes the domain extends from 6 — (pole) 
to — it/2 (equator). For odd I modes, the domain 
extends to 8 — ir. The number of grid points we employ 
is 200 x 80 for I being even and 160 x 120 for I being 
odd, in r and 6 respectively. The boundary conditions 
are implemented in the same way as in FSK. 

Our numerical evolution code was thoroughly 
tested in FSK, by comparing evolutions of perturbed 
spherical stars with results from perturbation theory 
obtained with an independent eigenvalue code. Since 
then, the code has been upgraded to run efficiently on 
a NEC SX-5/3C vector supercomputer. This modifica- 
tion was essential for doing a large number of numerical 
evolutions for many rotational periods and with a large 
number of grid-points. 



3 PERTURBATION OF THE INITIAL 
DATA 

The accurate computation of mode frequencies in a ro- 
tating star requires an appropriate excitation of the 
equilibrium initial data. When doing so, it is possible to 
obtain the frequencies of the excited modes with good 
accuracy, through a Fourier transform of the time evo- 
lution of the hydrodynamical variables, provided that 
the evolution time is much larger than the period of os- 
cillations. As in the hydrodynamical evolution we are 
using the 3+1 formulation (Banyuls et al 1997), the 
oscillation frequencies of the various evolved variables 
are obtained with respect to the coordinate time at a 
given location. This corresponds to the frequency of 
oscillations in a reference frame attached to an iner- 
tial observer at infinity. To increase the accuracy in 
the computation of the frequencies, we search for the 
zeros of the first derivative of the Fourier transform 
(with respect to the frequency) , using second-order ac- 
curate central differences. These zeros correspond to 
maxima in the Fourier transform, which (except for 
high-frequency noise) correspond to the excited modes 
of oscillation. This procedure is done at several points 
inside the star and the frequencies thus determined 
are found to be the same for each mode, i.e., all the 
modes that we identify are global discrete modes. For 
the resolution employed we estimate the accuracy of 
the computed frequencies to be of the order of 1% — 2%. 
The different pulsation frequencies are identified with 
specific normal modes by comparing frequencies of the 
non-rotating star to known eigenfrequencies from per- 
turbative normal-mode calculations. 
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Figure 1. 2-D greyscale plot of the initial data used for 
mode-excitation: density perturbation (I = 0) and vg (I ^ 
0). The darkest area corresponds to the maximum of the 
perturbation, while the lightest area corresponds to its min- 
imum values. The thick solid line indicates the location of 
the surface of the star. The depicted initial model is the 
fastest rotator of our sample in Table 1. 



As demonstrated in FSK, the small-amplitude pul- 
sations in the nonlinear, fixed spacetime evolutions cor- 
respond to linear normal modes of pulsation in the 
relativistic Cowling approximation (McDermott et al 
1983), in which perturbations of the spacetime are ig- 
nored. The existence of a numerical viscosity inher- 
ent to the numerical scheme, damps the pulsations of 
the star. Therefore, high-resolution grids are preferred 
to reduce the damping, especially for the higher fre- 
quency modes, which are damped faster. In addition, 
our numerical scheme requires the presence of a ten- 
uous, constant-density "atmosphere" surrounding the 
star, which is reset to its initial state after each time- 
step (in such a way that the stellar surface is allowed 
to contract or expand), introduces an additional nu- 
merical damping of the pulsations, due to the finite- 
differencing at the surface of the star. In order to min- 
imize this effect, the density of the "atmosphere" has 
to be small enough to be dynamically unimportant. As 
already mentioned, a typical value of 10~ 6 pc is appro- 
priate for this purpose. 

We use analytic eigenfunctions to excite particular 
oscillation modes. For the I = modes the initial equi- 
librium values of density and pressure, po and po, are 
perturbed to nonequilibrium values, p — po + Sp and 
p = po + Sp by the eigenfunctions 



where p c is the central density, r s (8) is the coordinate 
radius of the surface of the star (which depends on the 
polar angle 6) and T is the adiabatic index of the ideal 
gas EOS, p = (r — l)pe, related to the polytropic index 
by the equation V = 1 + 1/iV (for isentropic stars). The 
amplitude of the excitation, A, is typically chosen to 
be in the range 0.001 to 0.005. 

For the excitation of the I = 1, 2 and 3 modes we 
add a small non-zero 8- velocity component to perturb 
the initial vanishing value. More precisely, for I = 1 we 
have 



vg = A sin 



for I = 2 



vg = A sin 



rs(0) 



r.{9) 



and for I = 3 



vg — A sin ( —jtt ) sin 6(1 — 5 cos 2 6). 
r s (6) 



(5) 



(6) 



(7) 



Sp = Ap c sin 



(3) 



The particular radial dependence in the above eigen- 
functions is chosen so that the perturbations vanish 
at the surface of the rotating star. We have found that 
this is necessary for our numerical scheme to work. Fig- 
ure ^ displays a 2-D greyscale plot of the above pertur- 
bations, for the model near the mass-shedding limit. 
The darkest area corresponds to the maximum of the 
perturbation, while the lightest area corresponds to its 
minimum values. We note that all models have equa- 
torial plane symmetry with respect to the z = axis, 
even though we use a grid extending to 6 — tt for odd 
Z's. The choice of these trial eigenfunctions allows us to 
compute the four lowest frequencies quite accurately 
for all the considered I modes, as we discuss next. 



4 RESULTS 

We turn now to presenting our numerical results con- 
cerning the frequencies of axisymmetric pulsations of 
uniformly and rapidly rotating neutron stars. 



4.1 Quasi-radial modes 

Figure ^ shows the time evolution of the radial velocity 
component for the most rapidly rotating model in our 
sequence. The initial equilibrium model has been per- 
turbed with an / = perturbation according to Eqs. (^) 
and (^). The final evolution time corresponds to 10ms 
and the oscillations are measured at half the radius of 
the star and at an angle of 6 = 7r/4. The oscillatory 
pattern depicted in this figure is typical to all our sim- 
ulations: it is mainly a superposition of the lowest-order 
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Figure 2. Time evolution of the radial velocity of the most 
rapidly rotating model of our sequence, O = 0.522 X 10 4 s — 1 . 
An I = perturbation has been applied to the equilibrium 
data. The pulsations are mainly a superposition of the nor- 
mal modes of the star. 
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Figure 3. Logarithm of the amplitude (in arbitrary scale) 
resulting from the Fourier transform of the radial velocity 
evolution shown in Fig. ^ It is possible to identify in this 
plot the frequencies of the fundamental quasi-radial mode as 
well as up to three harmonics. Additionally, / and overtones 
of p— modes can also be identified. 



normal modes of the fluid. The amplitude of the oscil- 
lations is damped due to the inherent viscosity of the 
numerical scheme. The high frequency normal modes 
are usually damped faster and at the final time the star 
is pulsating mostly in a few lowest frequency modes. 

The frequencies of the axisymmetric modes are ob- 
tained by a Fast Fourier Transform of the time evo- 
lution of selected hydrodynamical variables (both the 
density and the components of the velocity). Figure H 



Table 2. Fundamental, first, second and third 
overtones (F, Hi, H2 and H3, respectively) of the 
quasi-radial (I = 0) modes for a sequence of ro- 
tating stars of same central density. The angu- 
lar velocity Si^ a t the mass-shedding limit is 
0.8094 X 10 4 s" 1 for this sequence. 



n 


F 


Hi 


H 2 


H 3 


(10 4 s-l) 


(kHz) 


(kHz) 


(kHz) 


(kHz) 


0.0 


2.706 


4.547 


6.320 


8.153 


0.218 


2.657 


4.467 


6.215 


8.005 


0.306 


2.619 


4.409 


6.202 


8.005 


0.371 


2.579 


4.385 


6.234 


8.096 


0.399 


2.553 


4.377 


6.243 


8.098 


0.423 


2.535 


4.371 


6.241 


8.134 


0.445 


2.510 


4.362 


6.266 


8.171 


0.465 


2.495 


4.356 


6.262 


8.171 


0.482 


2.476 


4.366 


6.274 


8.197 


0.498 


2.456 


4.357 


6.270 


8.130 


0.511 


2.442 


4.350 


6.297 


8.030 


0.522 


2.417 


4.337 


6.255 


7.987 



shows the Fourier transform of the radial velocity evo- 
lution depicted in Fig. |^. The main mode which is 
excited is the fundamental quasi-radial f-mode. Its 
higher harmonics (H1-H3) are also excited, as well as 
several other, non-radial modes. The amount of excita- 
tion of each mode depends on the correlation between 
the mode eigenfunction and the applied perturbation. 

As it is apparent in Figure ^, a dense spectrum 
of modes appears when one applies generic perturba- 
tions that do not correspond to the eigenfunctions of a 
particular normal mode only. Thus, in order to identify 
the peaks in the Fourier transform with specific normal 
modes, we rely on the previously known frequencies in 
the non-rotating limit (see FSK). As the rotation rate 
is increased, we follow the change in the location of the 
various peaks, keeping in mind that at large rotation 
rates apparent crossings of frequencies occur. In such a 
case, the amplitude of the Fourier transform at various 
points inside the star (which correlates with the mode 
eigenfunction) is used as a guide in deciding about the 
correct identification of the mode frequency. The spe- 
cific values for the frequencies of the fundamental F 
and higher harmonics Hi , H2 , -H3 quasi-radial modes 
for the sequence of rotating stars considered here, are 
shown in Table 0. 

As frequencies of quasi-radial modes have been 
computed previously by Yoshida & Eriguchi (2000), 
as an eigenvalue problem, we use those results to com- 
pare the values obtained with our code for a soft poly- 
trope, with p c = 8.1 x 10~ 4 , N = 1.5 and K = 
4.349. We have compared the models corresponding 
to = 0,7.1379 x 10~ 3 and 1.4094 x 10~ 2 , which, 
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Table 3. Comparison of I = quasi-radial pul- 
sation frequencies, obtained with the present non- 
linear evolution code, to linear perturbation mode 
frequencies in the relativistic Cowling approxima- 
tion (Yoshida & Eriguchi 2000). The equilibrium 
model is a N = 1.5, K = 4.349 relativistic poly- 
trope with p c = 8.1 X 10~ 4 . 





Mode 


Y&E 


present 


Difference 


(io- 3 ) 




(kHz) 


(kHz) 


(per cent) 





F 


1.674 


1.678 


0.3 




111 


2.758 


2.807 


1.7 




112 


3.793 


3.841 


1.3 


7.1379 


F 


1.646 


1.670 


1.5 




111 


2.696 


2.735 


1.4 




112 


3.728 


3.761 


0.9 


14.094 


F 


1.545 


1.553 


0.5 




111 


2.572 


2.595 


0.9 




112 


3.664 


3.642 


0.6 



in the notation of Yoshida & Eriguchi correspond to 
/rot = 0,0.1415 and 0.2794, respectively. The results 
of this comparison are presented in Table ^. We note 
that the agreement is very good, especially for the most 
rapidly rotating model, the differences always being be- 
low 2%. As there may be some small differences in the 
construction of the equilibrium model, the actual accu- 
racy of our code is better than the relative differences 
shown in Table til 



4.2 Non-radial modes 

In Figures (^j^) we plot the time evolutions of the po- 
lar velocity component, along with the corresponding 
Fourier transforms for non-radial I = 1,2 and 3 per- 
turbations, for the fastest rotating model in Table 1. 
The time evolutions shown are measured at half the 
star radius and 9 — 2tt/3 for I — 1,3, and 9 = 7r/3 
for 1 = 2. The above time evolutions show the same 
qualitative behavior already described for the 1 = 
modes. The lowest frequency (and dominant) I = 1 
mode we excite in our time-evolutions (labeled 1 /) has 
no nodes along the radial direction and behaves like a 
fundamental mode. We point out that the / = 1 fun- 
damental mode does not exist when one considers the 
full set of equations (i.e. including the perturbation of 
the metric), due to momentum conservation (it would 
correspond to a displacement of the center of mass). 
In the Cowling approximation, however, momentum is 
not conserved, as the perturbation in the metric is ne- 
glected. In this approximation, only the fluid oscillates 
in a fixed background metric and an oscillation of the 



center of mass is allowed, as the fixed metric acts as a 
restoring force. 

For I = 2, a larger number of modes is excited by 
the initial perturbation. At late times, the evolution 
is mainly a superposition of the fundamental 2 /-mode 
and the fundamental quasi-radial F-mode, as it is ev- 
ident from the amplitudes of the various modes in the 
corresponding Fourier transform. 

For I = 3, a large number of modes is also ex- 
cited by the initial perturbation, as in the I = 2 case. 
At late times, the evolution is mainly a superposition 
of the fundamental 3 f- and 1 /-modes and of the fun- 
damental quasi-radial F'-mode. The Fourier amplitude 
of the 1 /-mode is larger than the amplitude of the 
3 /-mode, which shows that the eigenfunction of the 
I = 3 modes near the mass-shedding limit is signifi- 
cantly modified by rotation, so that the / = 3 part of 
the 1 /-eigenfunction correlates better with the generic 
eigenfunction we used to excite I = 3 modes, than the 
3 /-eigenfunction itself. 

The frequencies of all identified / = 1,2 and 3 
modes are displayed in detail in Tables ^-^[ The modes 
are labeled as / for the fundamental modes and l p n for 
the p-modes of order n. A plot of all mode-frequencies 
as a function of the rotation rate is shown in Figure [n]. 
At rotation rates below ~ 50% of the mass-shedding 
limit, the frequencies of the lowest-order modes are 
not significantly affected by rotation. This is consis- 
tent with previous results in the slow-rotation approx- 
imation (Hartle & Friedman 1975) and in the New- 
tonian limit. For larger rotation rates, however, the 
I = and I = 1 overtones have a tendency to in- 
crease in frequency with rotation rate, while the I = 2 
and I = 3 overtones have the opposite tendency. As a 
result, several apparent mode crossings take place be- 
tween different modes. This has been observed before 
by Clement (1986) for Newtonian axisymmetric modes 
and by Yoshida & Eriguchi (2000) for the relativistic 
quasi-radial modes. 

In the above studies, when one follows an eigen- 
frequency continuously from the non-rotating limit to 
the large rotation rates, then, at apparent mode cross- 
ings, the continuous lines corresponding to different 
frequency sequences do not cross, which is normally 
called "avoided crossing". Such avoided crossings can 
occur in two ways: the eigenfunction along a continu- 
ous frequency sequence remains that of the same mode 
(which is the usual type of avoided crossing in non- 
rotating stars, see Unno et al 1989) or a different mode 
appears in the same continuous frequency sequence af- 
ter the avoided crossing. The latter case is encountered 
in rotating stars, as in Figure [Hj. To distinguish the 
two cases, we prefer to use the term "apparent cross- 
ing" for rotating stars (as was done in Clement 1986), 
which refers to a mode-sequence, rather than a fre- 
quency sequence. 

Finally, we note that in Clement (1981) it is 
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Figure 4. Time evolution of the polar velocity component 
for the fastest rotator, Q = 0.522 X 10 4 s _1 . An I = 1 per- 
turbation has been applied to the equilibrium data. 



Figure 6. Same as Fig. ^ but showing the time evolution 
for the I = 2 perturbation. 




4.0 6.0 
f(kHz) 



Figure 5. Fourier transform of the polar velocity evolution 
shown in Fig. |l| The different axisymmetric, £ = 1 modes 
are conveniently labelled. 



claimed that the axisymmetric 2 /-mode lies on a 
continuous branch, near the mass-shedding limit, for 
rapidly rotating Newtonian stars. In our relativistic 
computation, no such behavior was found for any of 
the modes studied. Within the numerical resolution 
employed all identified modes were found to be dis- 
crete. 



4.3 Inertial modes 

Apart from the quasi-radial and /- and p-modes, a 
number of axisymmetric inertial modes was also ex- 
cited in our numerical evolutions, which can be seen 



as low-frequency peaks in our Fourier transforms (Fig- 
ures |||9]). These modes exist in isentropic stars (such 
as those considered here) as a mixture of axial r-modes 
and polar (/-modes (see Lockitch & Friedman 1999). 
Nonaxisymmetric inertial modes have been computed 
as an eigenvalue problem for slowly-rotating relativis- 
tic stars (Lockitch 1999, Lockitch, Andersson & Fried- 
man 2000) , but frequencies for axisymmetric modes are 
not available yet. This makes their identification diffi- 
cult, as, in our simulations, many modes with similar 
frequencies appear. Furthermore, the spacing between 
these frequencies is of the same order as one would ex- 
pect the difference between the relativistic frequencies 
and the Newtonian frequencies (computed in Lockitch 
& Friedman 1999) to be. Therefore, it would be too 
venturous, at this point, to attempt an identification 
with specific normal modes, without prior knowledge 
of some of these frequencies in relativity, at least for 
slow-rotation. 



5 SUMMARY AND OUTLOOK 

We have presented a comprehensive study of all low- 
order axisymmetric modes of uniformly and rapidly 
rotating relativistic stars in the Cowling approxima- 
tion. This investigation has been carried out by numer- 
ically evolving initial perturbed equilibrium configura- 
tions with an axisymmetric, nonlinear, relativistic hy- 
drodynamics code. The simulations were performed us- 
ing a high-resolution shock-capturing finite-difference 
scheme accurate enough to maintain the initial rota- 
tion law for a sufficient number of rotational periods. 

Through Fourier transforms of the time evolution 
of selected fluid variables we computed the frequen- 
cies of non-radial, axisymmetric modes (with angular 
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Figure 7. Same as Fig. |] but showing the frequencies of Figure 9. Same as Fig. ^ but showing the frequencies of 

the 1 = 2 axisymmetric modes. the I = 3 axisymmetric modes. 
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Figure 8. Same as Fig. Q but showing the time evolution 
for the I = 3 perturbation. 



Table 4. /-mode and first three p-modes corre- 
sponding to I = 1. Further details in Tables [l] 
and El 



n 


7 


V 


V 


1 P3 


(10 4 s"l) 


(kHz) 


(kHz) 


(kHz) 


(kHz) 


0.0 


1.335 


3.473 


5.335 


7.136 


0.218 


1.349 


3.464 


5.318 


7.134 


0.306 


1.356 


3.453 


5.317 


7.152 


0.371 


1.364 


3.446 


5.320 


7.172 


0.399 


1.369 


3.442 


5.322 


7.193 


0.423 


1.371 


3.438 


5.325 


7.214 


0.445 


1.373 


3.434 


5.328 


7.238 


0.465 


1.375 


3.429 


5.333 


7.223 


0.482 


1.376 


3.422 


5.339 


7.349 


0.498 


1.376 


3.417 


5.340 


7.288 


0.511 


1.375 


3.407 


5.337 


7.281 


0.522 


1.375 


3.393 


5.335 


7.318 



momentum indices Z = 0, 1, 2 and 3) of rapidly rotat- 
ing stars. Therefore, we have extended previous results 
by Yoshida & Eriguchi (2000) which were mainly re- 
stricted to quasi-radial modes. We have presented re- 
sults for a complete sequence of rotating stars, ranging 
from the non-rotating case to rapid rotation near the 
mass-shedding limit. Apparent crossings between dif- 
ferent modes appear for rapidly rotating stars, as a 
result of the different influence of rotation on quasi- 
radial and Z = 1 modes than on modes with Z > 2. This 
different behavior may be related to the fact that the 
rotational deformation of the equilibrium star appears 
first as an Z = 2 term. Several axisymmetric inertial 
modes were also excited in our simulations. However, 
a definitive identification of the observed frequency 



peaks with specific modes will only be possible when 
mode-frequencies on the slow-rotation approximation 
be computed as an eigenvalue problem. Alternatively, 
a determination of mode-eigenfunctions in our simula- 
tions (and a comparison to mode-eigenfunctions in the 
Newtonian limit) may also allow the identification of 
such inertial modes. 

In following work we plan to study axisymmet- 
ric modes of differentially rotating stars with realistic 
equations of state. Moreover, the implementation of the 
hydrodynamic equations and numerical techniques em- 
ployed in the present work has been recently extended 
(Dimmelmeier et al 2000) to allow for gravitational 
field dynamics through the so-called conformally flat 



© RAS, MNRAS 000, 



Axisymmetric Modes of Rotating Relativistic Stars 9 



Table 5. /-mode and first three p-modes corre- 
sponding to I = 2. Further details in Tables []] 
and|[ 
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Table 6. /-mode and first three p-modes corre- 


sponding to / 


= 3. 


Further 


details in 


Tables |j 
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Figure 10. Frequencies of the lowest three quasi-radial modes vs. the ratio of angular velocitvof the star Q to the angular 
velocity at the mass-shedding limit Sljf, for the sequence of rotating relativistic stars in Table hi. 
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